Coalescing binaries as possible standard candles 
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Gravitational waves detected from well-localized inspiraling binaries would allow to determine, 
directly and independently, both binary luminosity and redshift. In this case, such systems could 
behave as "standard candles" providing an excellent probe of cosmic distances up to z < 0.1 and 
thus complementing other indicators of cosmological distance ladder. 
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I. INTRODUCTION 

A new type of standard candles, or, more appropri- 
ately, standard sirens, could be achieved by studying co- 
alescing binary systems [l], Q • These systems are usually 
considered strong emitter of gravitational waves (GW), 
ripples of space-time due to the presence of accelerated 
masses in analogy with the electromagnetic waves, due 
to accelerated charged. The coalescence of astrophysical 
systems containing relativistic objects as neutron stars 
(NS), white dwarves (WD) and black holes (BH) consti- 
tute very standard GW sources which could be extremely 
useful for cosmological distance ladder if physical features 
of GW emission are well determined. These binaries sys- 
tems, as the famous PSR 1913+16 0-01, have compo- 
nents that are gradually inspiralling one over the other 
as the result of energy and angular momentum loss due 
to (also) gravitational radiation. As a consequence the 
GW frequency is increasing and, if observed, could con- 
stitute a "signature" for the whole system dynamics. The 
coalescence of a compact binary system is usually clas- 
sified in three stages, which are not very well delimited 
one from another, namely the inspiral phase, the merger 
phase and the ring-down phase. The merger phase is 
the process that proceeds until the collision of the bod- 
ies and the formation of a unique object. Its duration 
depends on the characteristics of the originating stars 
and emission is characterized by a frequency damping 
in the time. In merger phase, stars are not modelled 
as rigid sphere due to the presence of a convulsive ex- 
change of matter ((J. GW emission from merger phase 
can only be evaluated using the full Einstein equations. 
Because of the extreme strong field nature of this phase, 
neither a straightforward application of post-Newtonian 
theory nor any perturbation theory is very useful. Re- 
cent numerical work has given some insight into 
the merger problem, but there are no reliable models for 
the waveform of the merger phase up to now. Gravita- 
tional radiation from the ring-down phase is well known 
and it can be described by quasi- normal mode s II Oil . The 
relevance of ring-down phase is described in [ll| . Tem- 
poral interval between the inspiral phase and the merger 



one is called coalescing time, interesting for detectors as 
the American LIGO (Laser Interferometer Gravitational- 
Wave Observatory) and French/Italian VIRGO [lj]. 
Coalescence is a rare event and therefore to see several 
events per year, LIGO and/or VIRGO must look far be- 
yond our Galaxy. For example, the expected rate for 
a NS-NS system (determined from the observed popu- 
lation of NS-NS binaries) is found to be 80±7q° Myr" 1 
per galaxy [14[. From this figure, one finds that the ex- 
pected rate for the today available sensitivities of today 
LIGO and VIRGO is of the order 35±|g x 1(T 3 yr" 1 , 
while for the advanced version of such interferometers, 
the rate is more interesting being 190^150 y r_1 j that is 
the probability ranges from one event per week to two 
events per day. A remarkable fact about binary coales- 
cence is that it can provide an absolute measurement of 
the source distance: this is an extremely important event 
in Astronomy. In fact, for these systems, the distance 
is given by the measure of the GW polarization emit- 
ted during the coalescence. One of the problems which 
affects the utilization of coalescing binaries as standard 
candles is the measure of the redshift of the source as well 
as the measure of the GW polarization (at present GWs 
have not been still experimentally observed). A solution 
for the redshift determination could be the detection of 
an electromagnetic counterpart of the coalescing system 
(e.g. the detection of an associated gamma ray burst) 
or the redshift measurement of the host galaxy or galaxy 
cluster at their barycenter [ll| . Recent evidence supports 
the hypothesis that many short-hard gamma-ray bursts 
could be associated with coalescing binary systems in- 
deed (l6l - fli| . In this paper, we want to show that such 
systems could be used as reliable standard candles. 

In Sect. HQ we briefly sketch the GW emission from 
coalescing binary systems in circular orbit at cosmologi- 
cal distance. In Sect. IIII1 we simulate various coalescing 
binary systems (WD- WD, NS-NS, BH-BH) at redshift 
z < 0.1 because of the observational limits of ground- 
based-interferometers as LIGO and VIRGO. The goal of 
the simulation is the measure of the Hubble constant and 
consequently the use of these systems as standard can- 
dles [20] - This new type of standard candles will be able 
to increase the confidence level on the other " traditional" 
standard candles in Astronomy and, moreover, it could 
constitute an effective tool to measure distances at larger 
redshifts. Conclusions are drawn in Sect. IIVI 
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II. GRAVITATIONAL RADIATION FROM A 
COALESCING BINARY SYSTEM 

For a detailed exposition of GW theory see, e.g. [2lT — 
[23j . Let us consider here an isolated, far away and slowly 
moving source. In this approximation, we can write the 
GW solution as: 



2G d 2 I k 
r dr' 2 



(1) 



where is the quadrupole momentum tensor of the en- 
ergy density of source, conventionally defined as: 



Iij(t)= / yVT 00 (t,y)d 3 y, 



(2) 



a tensor defined at any constant time surface, and t R is 
the retarded time [2l[ . The distance between source and 
observer is denoted as r. 

The gravitational wave produced by an isolated non- 
relativistic object is therefore proportional to the second 
derivative of the quadrupole momentum of the energy 
density at the point where the past light cone of the 
observer intersects the source. In contrast, the leading 
contribution to electromagnetic radiation comes from the 
changing dipole momentum of the charge density. 

The above result can be specified for two point masses 
M\ and M2 in a circular orbit. In the quadrupole ap- 
proximation, the two polarization amplitudes of GWs at 
a distance r from the source are given evaluating Eq.([T]) 
to lowest order in v/c, that is: 




Figure 1: Lines of force associated to the + (left panel) and 
x (right panel) polarizations. 



$ is given by the expression: 



$(t) = 2tt J dt'f{t') . 



(6) 



'0 



Note that for a fixed distance r and a given frequency 
/, the GW amplitudes are fully determined by /jM 2 ' 3 = 
Mc 5 ^ 3 , where the combination: 



Ma = 



(MiM 2 ) 3/5 
(Mi +M 2 ) 1/5 ' 



(7) 



Mo 



is called chirp mass of the binary, here M = Mi 

MiM 2 

is the total mass of the system and u = is the 

p Mi+M 2 

reduced mass of the system. Introducing the coalescence 



time t = t c 



t and integrating Eq. ([5]), we get: 



h+ (t) 



4 /&?./< 



G 



5/3 



2/3 



cos [$ {t R )\ 



(3) 



h x (t) 



= 4 / GMc \ 
r \ c 2 / 



5/3 



7T/ (ti 



2/3 



cosz 



sin [$ (t R )} 



(4) 

where i is the binary inclination angle such that i = 
90° corresponds to a system visible edge-on. These are 
traditionally labeled " plus" and " cross" from the lines of 
force associated with their tidal stretch and squeeze (see 
Fig. [T]). Here / is the frequency of the emitted GWs 
(twice the orbital frequency) . The rate of the frequency 
change is 22] : 



/ = £\b/3 (GMc\ 5/3 fll/3 



(5) 



where t R is the so called "retarded time" and the phase 



/~130 



1.21M 



M, 



c 



5/8 



1 sec 



3/8 



Hz. 



(8) 



Eq. (jHJl predicts coalescence times of t ~ 
17mm, 2sec, 1msec for / - 10, 100 , lOOOiiz. 

After averaging over the orbital period and the ori- 
entations of the binary orbital plane, one arrives at the 
average (characteristic) GW amplitude: 



h(f,M c ,r) = ((hl) + (hl)) 1/2 = 



32 \ 1/2 G 5 ^ 3 M C /3 



(*/) 



2/3 



(9) 



For a binary at the cosmological distance, i.e. at red- 
shift z where GWs propagate in a Friedmann- Robertson- 
Walker Universe, these equations arc modified in a very 
straightforward way: 

• The frequency that appears in the above formulae 
is the frequency measured by the observer, fobs, 
which is red-shifted with respect to the source fre- 
quency f s , i.e. f obs = f s /(l + z), and similarly t 
and t R are measured with the observer clocks. 
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• The chirp mass Mc has to be replaced by Aic = 
M c {l + z). 

• The distance r to the source has to be replaced by 
the luminosity distance cLl{z). 

Inserting the following quantity: 



hAt) *(^Y 3 (im)'<\ (10) 



d L (z) V c 2 



we can rewrite the expressions for the polarization "+" 
and " x" as: 



1 -I- f*ns 7 

h + (t) = h c (t R )—- cos [$(*«)], (11) 



and 



h x {t) = h c (t fl )cos«sin[$ (t R )} 



(12) 



Explicating the dependence on the chirp mass redshift, 
we can obtain the luminosity distance using the equa- 
tion (p~0|) linked directly to the GW polarization: 



d L (z) 



h c (t) 

4(1 + zf/ 3 
h c (t) 



GM c (l + z) 



5/3 



[*/(*)] 



2/3 



GM< 



5/3 



(13) 



[*/(*)] 



2/3 



where H(z) is the value of the Hubble parameter at red- 
shift z. Knowing di,{z), we can therefore obtain H(z). 
This shows that the luminosity distance function dziz) 
is an extremely important quantity, which encodes the 
whole expansion history of the Universe. Coalescing bi- 
naries could be standard candles (or precisely standard 
sirens) in the following sense. Suppose that we can mea- 
sure the amplitudes of both polarizations h + , h x , as well 
as fobs- From the ratio of h + and h x , we can obtain 
the value of the inclination of the orbit, besides, evalu- 
ating f b s at a given frequency, we can obtain A4c- If 
we are capable of measuring the redshift z of the source, 
we have found a gravitational standard candle since we 
can obtain the luminosity distance from Eq. (JT3J) and 
then evaluate the Hubble constant Hq. The difference 
between gravitational standard candles and the "tradi- 
tional" standard candles is that the luminosity distance 
is directly linked to the GW polarization and there is 
no theoretical uncertainty on its determination a part 
the redshift evaluation. Various possibilities have been 
proposed. Among these there is the possibility to see 
an optical counterpart. In fact, it can be shown that 
observations of the GWs emitted by inspiralling binary 
compact systems can be a powerful probe at cosmological 
scales. In particular, short GRBs appear related to such 
systems and quite promising as potential GW standard 
sirens [34]). On the other hand, the redshift of the binary 
system can be associated to the barycenter of the host 
galaxy or the galaxy cluster as we are going to do here. 



III. NUMERICAL SIMULATION 



Let us recall that the luminosity distance dj, of a source 
is defined by 



(14) 



where J- is the flux (energy per unit time per unit area) 
measured by the observer, and C is the absolute luminos- 
ity of the source, i.e. the power that it radiates in its rest 
frame. For small redshifts, dL is related to the present 
value of the Hubble parameter Hq and to the deceleration 
parameter go by 



H od L 1 2 

=z + -(l- q )z + 

c 2 



(15) 



The first term of this expansion gives the Hubble law 
z ~ (_ffo/c)di, which states that redshift is proportional 
to the distance. The term 0(z 2 ) is the correction to the 
linear law for moderate redshifts. For large redshifts, 
the Taylor series is no longer appropriate, and the whole 
expansion history of the Universe is encoded in a function 
di(z). For a spatially flat Universe, one finds 



d L {z) = c(l + z) 



dz' 



(16) 



We have simulated various coalescing binary systems 
at redshifts z < 0.1 In this analysis, we do not consider 
systematic errors and errors on redshifts to obviate the 
absence of a complete catalogue of such systems. The 
choice of low redshifts is due to the observational limits 
of ground-based interferometers like VIRGO or LIGO. 
Some improvements are achieved, if we take into ac- 
count the future generation of these interferometers as 
Advanced VIRGO (35|and Advanced LIGO 36]. Ad- 
vanced VIRGO is a major upgrade, with the goal of 
increasing the sensitivity by about one order of mag- 
nitude with respect to VIRGO in the whole detection 
band. Such a detector, with Advanced LIGO, is ex- 
pected to see many events every year (from 10s to 100s 
events/year). For example, a NS-NS coalescence will be 
detectable as far as 300 Mpc. In the simulation pre- 
sented here, sources are slightly out of LIGO- VIRGO 
band but observable, in principle, with future interfer- 
ometers. Here, we have used the redshifts taken by 
NASA/IPAC EXTRAGALACTIC DATABASE H3, and 
we have fixed the redshift using z at the barycenter of the 
host galaxy/cluster), and the binary chirp mass Mq, typ- 
ically measured, from the Newtonian part of the signal 
at upward frequency sweep, to ~ 0.04% for a NS/NS bi- 
nary and ~ 0.3% for a system containing at least one BH 
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[28l |38||. The distance to the binary ("luminosity dis- 
tance" at cosmological distances) can be inferred, from 
the observed waveforms, to a precision ~ 3/p < 30%, 
where p = S/N is the amplitude signal-to-noise ratio in 
the total LIGO network (which must exceed about 8 in 
order that the false alarm rate be less than the threshold 
for detection). In this way, we have fixed the character- 
istic amplitude of GWs, and frequencies are tuned in a 
range compatible with such a fixed amplitude, then the 
error on distance luminosity is calculated by the error on 
the chirp mass with standard error propagation. 

The systems considered are NS-NS, BH-BH and WD- 
WD. For each of them, a particular frequency range and 
a characteristic amplitude (beside the chirp mass) are 
fixed. We start with the analysis of NS-NS systems 
[Mc = 1-22M ) with characteristic amplitude fixed to 
the value 10~ 22 . In Table HI we report the redshift, the 
value of he and the frequency range of systems analyzed. 



Object 


z 


h c 


Freq. (Hz) 


NGC 5128 


0.0011 


1Q -22 


0- 


- 10 


NGC 1023 Group 


0.0015 


1Q -22 


0- 


- 10 


NGC 2997 


0.0018 


1Q -22 


5- 


- 15 


NGC 5457 


0.0019 


10~ 22 


10^20 


NGC 5033 


0.0037 


1Q -22 


25^35 


Virgo Cluster 


0.0042 


1Q -22 


30-^40 


Fornax Cluster 


0.0044 


1Q -22 


35-45 


NGC 7582 


0.0050 


10 -22 


45^55 


Ursa Major Groups 


0.0057 


1Q -22 


50^60 


Eridanus Cluster 


0.0066 


10~ 22 


55^65 



Table I: Redshifts, characteristic amplitudes, frequency range 
for NS-NS systems. 

In Fig. [21 the derived Hubble relation is reported, and 
in Fig. |3] residuals are reported. 

The Hubble constant value is 72 ± 1 km/sMpc in 
agreement with the recent WMAP estimation (Wilkin- 
son Microwave Anisotropy Probe [25]). The same proce- 
dure is adopted for WD- WD systems (M c = 0.69M Q , 
h c = 1CT 23 ) and BH-BH systems (M c = 8.67M , 
h c = 1CT 21 ). In Tables |H] and Hill we report the redshift, 
the value of he and the frequency range for BH-BH and 
WD- WD systems respectively. These simulations are re- 
ported in Fig. 2] and in Fig. [5] the corresponding residu- 
als have been checked and their goodness is the same of 
previous case so we omit them. 

For these simulations, the Hubble constant value is 
69 ± 2 km/sMpc and 70 ± 1 km/sMpc for BH-BH 
and WD- WD systems respectively, also in these cases 
in agreement with WMAP estimation. 



IV. CONCLUSIONS 

In this paper, we have considered simulated binary sys- 
tems whose redshifts can be estimated considering the 
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Figure 2: Luminosity distance vs redshift for simulated NS- 
NS systems. 
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Figure 3: Residuals of the plot luminosity distance vs redshift 
for simulated NS-NS systems. 



barycenter of the host astrophysical system as galaxy, 
group of galaxies or cluster of galaxies. In such a way, 
the standard methods adopted to evaluate the cosmic 
distances (e.g. Tully-Fisher or Faber- Jackson relations) 
can be considered as "priors" to fit the Hubble rela- 
tion. We have simulated various situations assuming 
NS-NS, BH-BH, and WD- WD binary systems. Clearly, 
the leading parameter is the chirp mass Mc, or its red- 
shifted counter-part Aic, which is directly related to the 
GW amplitude. The adopted redshifts are in a well- 
tested range of scales and the Hubble constant value is in 
good agreement with WMAP estimation. The Hubble- 
luminosity-distance diagrams of the above simulations 
show the possibility to use the coalescing binary systems 
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Object 


z 


h c 


Freq. (Hz) 


Pavo-Indus Sup. Cluster 


0.015 


io- 21 


65^70 


Abell 569 Cluster 


0.019 


io- 21 


75^-80 


Coma Cluster 


0.023 


io- 21 


100 -r 105 


Abell 634 Cluster 


0.025 


IO" 21 


110-M15 


Ophiuchus Cluster 


0.028 


io- 21 


130 ~ 135 


Columba Cluster 


0.034 


io- 21 


200 ^ 205 


Hercules Sup. Cluster 


0.037 


io- 21 


205 210 


Sculptor Sup. Cluster 


0.054 


ht 21 


340 ~ 345 


Pisces-Cetus Sup. Cluster 


0.063 


io- 21 


420 425 


Horologium Sup. Cluster 


0.067 


IO" 21 


450 ^ 455 



Object 


z 


h c 


Freq. (Hz) 


Eridanus Cluster 


0.0066 


1Q -23 


5-rl0 


Hydra Cluster 


0.010 


1Q -23 


15- 


-20 


Payo- Indus Sup. Cluster 


0.015 


1Q -23 


35- 


-40 


Perseus- Pisces Sup. Cluster 


0.017 


IO -23 


40- 


-45 


Abell 569 Cluster 


0.019 


1Q -23 


45- 


-50 


Centaurus Cluster 


0.020 


1Q -23 


45- 


-50 


Coma Cluster 


0.023 


1Q -23 


55- 


-60 


Abell 634 Cluster 


0.025 


1Q -23 


60- 


-65 


Leo Sup. Cluster 


0.032 


1Q -23 


85- 


-90 


Hercules Sup. Cluster 


0.037 


1Q -23 


100- 


- 105 



Table II: Redshifts, characteristic amplitudes, frequency 
range for BH-BH systems. 



Table III: Redshifts, characteristic amplitudes, frequency 
range for WD- WD systems. 
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Figure 4: Luminosity distance vs redshift for simulated BH- 
BH systems. 



Figure 5: Luminosity distance vs redshift for simulated WD- 
WD systems. 



as distance indicators and, possibly, as standard candles. 
The limits of the method are, essentially, the measure 
of GW polarizations and redshifts. Besides, in order to 
improve the approach, a suitable catalogue of observed 
coalescing binary-systems is needed. This is the main dif- 
ficulty of the method since, being the coalescence a tran- 
sient phenomenon, it is very hard to detect and analyze 
the luminosity curves of these systems. Furthermore, a 
few simulated sources are out of the LIGO- VIRGO band. 

Next generation of interferometer (as LISA (2(| or 
Advanced- VIRGO and LIGO) could play a decisive role 
to detect GWs from these systems. At the advanced 
level, one expects to detect at least tens NS-NS coalesc- 
ing events per year, up to distances of order 2 Gpc, mea- 
suring the chirp mass with a precision better than 0.1%. 
The masses of NSs are typically of order 1.4M©. Stellar- 
mass BHs, as observed in X-ray binaries, are in general 
more massive, typically with masses of order 10M©, and 
therefore are expected to emit even more powerful GW 



signals during their inspiralling and coalescing phases. 
The coalescence of two BHs, each one with 10 M©, could 
be seen by Advanced- VIRGO and Advanced-LIGO up to 
redshifts z ~ 2 — 3. [28] Furthermore, the LISA space in- 
terferometer, which is expected to fly in about 10 years, 
will be sensitive to GWs in the mHz region, which cor- 
responds to the wave emitted by supermassive BHs with 
masses up to 10 6 M©. Nowadays, supermassive BHs with 
masses between 10 6 and 10 9 M© are known to exist at 
the center of most (and probably all) galaxies, including 
our Galaxy. The coalescence of two supermassive BHs, 
which could take place, for instance during the collision 
and merging of two galaxies or in pre-galactic structure 
at high redshifts, would be among the most luminous 
events in the Universe. Even if the merger rate is poorly 
understood, observations from the Hubble Space Tele- 
scope and from X-ray satellites such as Chandra [13] have 
revealed that these merging events could be detectable 
at cosmological distances. LISA could detect them up 
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to z »j 10, [2^, H3| and it is expected to measure sev- 
eral events of this kind. The most important issue that 
can be addressed with a measure of oIl(z) is to under- 
stand "dark energy" , the quite mysterious component of 
the energy budget of the Universe that manifests itself 
through an acceleration of the expansion of the Universe 
at high redshift. This has been observed, at z < 1.7, us- 
ing Type la supernovae as standard candles [3l|, HH . A 
possible concern in these determinations is the absence 
of a solid theoretical understanding of the source. After 
all, supernovae are complicated phenomena. In partic- 
ular, one can be concerned about the possibility of an 
evolution of the supernovae brightness with redshift, and 
of interstellar extinction in the host galaxy leading to 
unknown systematics. GW standard candles could lead 
to completely independent determinations, and comple- 



ment and increase the confidence of other standard can- 
dles, (33j |. as well as extending the result to higher red- 
shifts. In the future, the problem of the redshift could 
be obviate finding an electromagnetic counterpart to the 
coalescence and short GRBs could play this role. 

In summary, this new type of cosmic distance indi- 
cators could be considered complementary to the tra- 
ditional standard candles opening the doors to a self- 
consistent gravitational astronomy. 
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